Dark matter and pulsar signals for 
Fermi LAT, PAMELA, ATIC, HESS and WMAP data 



V. Barger 1 , Y. Gao 1 , W.-Y. Keung 2 , D. Marfatia 3 , G. Shaughnessy 4 ' 5 
1 Department of Physics, University of Wisconsin, Madison, WI 53706 
2 Department of Physics, University of Illinois, Chicago, IL 60607 
3 Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045 
4 Department of Physics and Astronomy, Northwestern University, Evanston, IL 60208 and 
5 High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439 

We analyze new diffuse gamma-ray data from the Fermi Gamma-ray Space Telescope, which do not 
confirm an excess in the EGRET data at galactic mid-latitudes, in combination with measurements 
of electron and positron fluxes from PAMELA, Fermi and HESS within the context of three possible 
sources: dark matter (DM) annihilation or decay into charged leptons, and a continuum distribution 
of pulsars. We allow for variations in the backgrounds, consider several DM halo profiles, and account 
for systematic uncertainties in data where possible. We find that all three scenarios represent the 
data well. The pulsar description holds for a wide range of injection energy spectra. We compare 
with ATIC data and the WMAP haze where appropriate, but do not fit these data since the former 
are discrepant with Fermi data and the latter are subject to large systematic uncertainties. We 
show that for cusped halo profiles, Fermi could observe a spectacular gamma-ray signal of DM 
annihilation from the galactic center while seeing no excess at mid-latitudes. 

Introduction. Recent observations of an excess in cosmic ray positrons have led to much excitement in the particle 
physics community. At the heart of this excitement is the possibility that these may be a signal of the annihilation 
or decay of cold dark matter (DM) - the nonbaryonic, nonluminous component of matter that constitutes 20% of the 
energy budget of the universe, and is either exactly stable or quasi-stable over a time-scale that is much longer than 
the age of the Universe. 

The Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) experiment [ij reports 
a prominent upturn in the positron fraction from 10 — 100 GcV, in contrast to what is expected from high-energy 
cosmic rays interacting with the interstellar medium. (We employ the updated PAMELA spectrum derived 
using 2.5 times more data and refined background rejection methods, and which is slightly softer than the published 
spectrum in Ref. [l[0 This result confirms excesses seen in previous experiments, such as High Energy Antimatter 
Telescope (HEAT) [|| and Alpha Magnetic Spectrometer (AMS-01) 0]. The PAMELA satellite has a calorimeter 
and magnetic spectrometer calibrated up to 200 GeV and requires discrimination between positrons from protons at 
the level of 2 x 10~ 4 . Due to its small size, a large fraction of the electromagnetic shower from e^- is not contained, 
making the separation from protons and antiprotons (which do not shower) difficult. If the data are taken at face 
value, the rise in e + fraction is indicative of some source that injects high energy positrons that are not included 
in the standard astrophysical background estimates. PAMELA did not find an excess of antiprotons which may 
be indicative of a le ptop hilic DM particle [|| or pulsars 0, H, Q ■ Data obtained by the Advanced Thin Ionization 
Calorimeter (ATIC) [10( and the Polar Patrol Balloon and Balloon borne Electron Telescope with Scintillating fibers 
(PPB-BETS) [ll| experiments show an excess in the e + + e~ flux from 200 to 800 GeV. The challenge for these 
experiments, which can not discriminate charge, is to distinguish e^- from neutral pions at the level of 10~ 3 . The 
excess seen by these experiments is in an uncalibrated energy regime above 200 GeV. It is worth mentioning that ATIC 
has also found bumps in the spectra of light cosmic nuclei that remain puzzling [12]. Preliminary results from ATIC- 
4, which has a larger calorimeter to contain the electromagnetic showers, confirm their previous data. In addition, 
the High Energy Stereoscopic System (HESS) 7— ray telescope has inferred a flat but statistically limited e + + e~ 
spectrum between 340 GeV and 1 TeV [H[ which falls steeply above 1 TeV [l4| indicating a high energy limit of the 
new e^- source. The extraction of the e spectrum, which is the first of its kind, is made difficult by contamination 
from known 7— ray sources and the hadronic cosmic ray background. Although designed as a 7— ray observatory, 
the Large Area Telescope (LAT) of the Fermi Gamma-ray Space Telescope [15| can also detect e ± without charge 
discrimination. Recent high-statistics data from Fermi show an almost flat e + + e~ spectrum between 20 GeV and 
1 TeV, with a slight excess between 200 GeV and 1 TeV, but no bump- like structure jig , in contradiction with ATIC. 
It has been noted that a slightly harder supernova injection spectrum for the electron background can reproduce the 
Fermi data (at the expense of agreement with lower energy data), but can not account for the PAMELA excess [l7l |. 

Besides direct observations of excessive cosmic e , the Energetic Gamma Ray Experiment Telescope (EGRET) 
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observations provide hints of an excess of diffuse 7— rays from 10 — 50 GeV. Reanalysis of the data using improved 
sensitivity estimates [l8[ found a harder spectrum at these energies thanpreviously derived. The Fermi LAT 
observes 7— rays from 20 MeV to 300 GeV. Preliminary results from Fermi [20( show no readily discernible excess 
above backgrounds in the GeV range. 

The Wilkinson Microwave Anisotropy Probe (WMAP) observes microwave radiation in the 23 to 94 GHz range. 
A residual microwave background has been identified in WMAP data after known astrophysical backgrounds are 
subtracted from the skymap [2ll . [22j ■ This haze seems to originate from the synchrotron radiation of high energy 
electrons that may be attributable to annihilations of a thermal relic particle [23| . However, the systematic uncert ainty 
of the WMAP haze has not been robustly established, especially in the inner galactic region. As pointed out in Ref. [22| . 




the significance of the haze can be substantially reduced by allowing a spatial variation in the frequency dependence 
of the synchrotron radiation in the inner and outer galactic regions. 

In this letter, we perform a comparative study of three scenarios (described in the next section) that have been 
proposed to explain the aforementioned data. We do not include ATIC data in our statistical analyses since they are 
not corroborated by Fermi and we do not include the low energy HESS data sample [HI since the Fermi dataset is 
significantly larger in the overlapping energy region below 1 TeV. Nevertheless, we comment on ATIC data throughout 
the text. A distinguishing feature of our analysis is that we vary the power-law index for the electron background 
injection spectrum, account for cosmic ray propagation uncertainties and energy calibration uncertainties in data 
(which can rigidly shift spectra significantly), and consider several DM halo profiles. All these have the potential to 
alter conclusions drawn from less thorough analyses. 

Sources, (i) Annihilating Dark Matter. Genuinely stable DM, with a discrete symmetry that enforces stability, 
is central to many scenarios such as Supersymmetry, minimal Universal Extra Dimensions and Little Higgs with T- 
parity. Typically, the DM candidate interacts weakly and has a mass Mom in the 100 — 1000 GeV range. Throughout, 
we set the nonrelativistic thermally averaged annihilation cross section to the canonical value (av) = 3 x 10 -26 
cm 3 s _1 , required for saturating the DM relic density in the standard cosmological scenario. Direct annihilation or 
fragmentation to antiparticle states or high energy gamma-rays may yield a contribution above the astrophysical 
background. For recent studies related to the PAMELA and ATIC data see Refs. [241 [25l [26j. 

(ii) Decaying Dark Matter. If the exact discrete symmetry that protects the stability of the DM particle is slightly 
broken, DM may decay on time-scales longer than the age of the universe. Phenomenologically, this case is similar 
to DM annihilation with the exception that the DM mass is a factor of two larger, and the injection rate depends 
linearly (instead of quadratically) on the density of DM at the source. See Ref. [27]] for recent work. 

(in) Pulsars. An alternative explanation of the PAMELA and ATIC excess is pulsars that emit high energy 
positrons 7, 8, 9ij . Electrons in the intense rotating magnetic field that surrounds the neutron star can emit synchrotron 
radiation that is energetic enough to produce electron and positron pairs. The resulting positrons can propagate away 
from the pulsar via open magnetic field lines or via the pulsar wind. The resulting positron spectra can be modeled 
as a product of a power law and a decaying exponential with a cutoff at energy E p . 

To conveniently present our results in terms of their dependence on a generic energy scale, E s , of injected positrons 
for the three cases above, we define 



Cosmic ray species. Numerous astrophysical sources produce high energy cosmic rays that propagate through 
turbulent galactic magnetic fields. We utilize the GALPROP code [28| to simulate this propagation. For the distribu- 
tion of the DM halo we consider the isothermal (29|, Navarro, Frenk & White (NFW) [30| , Einasto [3l[ and Moore [32j 
profiles; see Fig. [T] 

Antimatter. Of the various cosmic ray species, antimatter provides one of the best avenues for probing DM because 
of relatively low backgrounds. The majority of the galactic e + background is created by the decay of cosmic rays or 
their interaction with the interstellar medium. 

(i) Annihilating Dark Matter. The positron spectrum from DM annihilation strongly depends on the underlying 
model and spin of the DM particle. To explain the hard positron spectra observed by PAMELA, a dominant annihila- 
tion fraction to charged leptons is preferred [24],[25|]. For the excess in the ATIC data, Mdm between 600 — 1000 GeV 
is required. However, the cross section needed to reproduce the relic density is about two orders of magnitude lower 
than that required to fit the e ± data. This has led many to invoke a boost factor from s-channel resonances [33j . 
the nonrelativistic enhancement from the Sommerfeld effect [34[, or DM clumping in the halo substructure [Hj]. A 
boost factor of astrophysical origin is easily motivated since local clumps in the halo distribution can enhance the flux 
from DM annihilation. N-body simulations suggest boost factors of order 10 [36], but it must be borne in mind that 
the simulations neglect ordinary matter. In principle, the flux from e ± , antiprotons and photons may have different 





Annihilating DM 
Decaying DM 
Pulsars 



(1) 



3 



10.0 



7.0 



5.0 



S 

u 

> 

O 3.0 



2.0 



1.5 



1.0 








, " 










-*:>.>> 






Einasto 




NFW 




_ Isothermal 





0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 
r [kpc] 



FIG. 1: rp(r) for the Moore, Einasto, NFW, and isothermal profiles with density 0.3 GeV/cm at 8.5 kpc. The analytical 
expressions are provided in the appendix. 



boost factors that are energy dependent, and depend on astrophysical assumptions such as the inner density profiles 
of clumps. We assume equal energy-independent boost factors for all species. 

In contrast to their positron data, PAMELA observes no excess in the antiproton spectrum and therefore rules 
out most processes that produce e + and hardonic states at roughly the same rate. For instance W, Z an d q uark 
pair productions generally produce too many antiprotons to agree with both the e + and p PAMELA data [24 |25| . 
However, if the DM particles annihilate dominantly into leptons, consistency with the antiproton data can be achieved 
even if the DM mass is a few hundred GeV due to the suppressed fragmentation to p. We consider cases in which 
DM annihilates to e + e~, and t + t~ with 100% branching fractions, and to these channels with equal branching 

fractions. 

(ii) Decaying Dark Matter. For fermion DM, decay channels such as W /± £ T with I a lepton, are possible, and the 
e + spectrum originates from equally weighted £ + and W fragmentation. However, the associated may contribute 
too many antiprotons to be consistent with PAMELA. Also, the 3-body decay mode t + l~v is permissible, but model- 
dependent. It is noteworthy that decay modes that are helicity suppressed are not intrinsically disfavored since the 
lifetime suggested by the data is significantly long. For the sake of direct comparison with DM annihilation, we restrict 
our study to 2-body decay modes into lepton pairs with equal branching fractions. Thus, for DM decay we effectively 
assume that the DM particle is a boson. 

A nice feature of DM decay is that no boost factors are required since the DM relic abundance is produced via 
annihilation while the astrophysical signals arise from DM decay. 

(Hi) Pulsars. It is unlikely that a single pulsar is responsible for both the PAMELA and ATIC results. The 
PAMELA data favor a mature pulsar (2 Myr) at a distance of a kpc, or a younger and closer pulsar like Geminga 
(0.37 Myr old at 160 pc), while the ATIC data require a very powerful and younger (0.5 Myr) pulsar at a distance of 
1 — 2 kpc Q. Rather than considering known nearby pulsars, we follow Ref. [37| and model the pulsar sources as a 
continuum distribution of pulsars throughout the galaxy in cylindrical coordinates as 



b(r-r ) 

p(r) = N [ — ) < 



(2) 



where zq = 0.2 kpc, r is the cylindrical radius from the galactic center parallel to the galactic disk, and r Q = 8.5 kpc. 
For mature pulsars, we take a = 1.0 and b — 1.8 [37]. N is a normalization that is determined by fitting PAMELA, 
Fermi and HESS data. 

We take the positron injection spectrum to be 
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where we fix a = 1.5 and allow E p to span the range from 150 GeV to 2 TeV. By fixing a, the number of indepen- 
dent parameters that describe the combined signal and background spectra for pulsars is the same as that for DM 
annihilations and decays. Note that variations in a can be somewhat compensated by N, thus making our results 
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only mildly sensitive to a over acceptable ranges. However, it should be noted that pulsars of various ages may well 
possess a continuum of cut-off energies Q and an age distribution may lead to nontrivial structures in the high energy 
spectrum, with young pulsars contributing more dominantly to the high energy range than mature pulsars. 

Pulsars do not produce antiprotons because of the absence of hardonic showers, and thus naturally explain why the 
PAMELA antiproton spectrum matches the expected spectrum. 

Light nuclei. Although not as compelling as the positron excess, the ATIC experiment has found excesses in the 
spectra of several light nuclei including C, N, O and Si, in the same energy range as the e + + e~ excess [12j. None of 
the scenarios under consideration offer an explanation of these anomalies, and we do not consider them any further. 

Photons. Under our assumption that the DM and pulsar sources produce only charged leptons directly, the main 
contribution to the flux of high energy photons is final state radiation (FSR) from the e + e~ mode and secondary 
decays. In addition, there is a continuum at lower energies from bremsstrahlung as the leptons leave the source. Once 
the leptons are ejected into the galactic environment, they interact with the galactic magnetic field and the interstellar 
radiation field (ISRF). They lose energy via inverse compton (IC) scattering which produces gamma radiation and 
via synchroton radiation as radiowaves. The diffuse gamma-ray background can be modeled as a sum of galactic 
pion decay (which dominates in the 1 — 10 GeV range), IC scattering, bremsstrahlung and the isotropic extragalactic 
(EG) HU diffuse gamma-rays (estimated up to about 50 GeV and which become increasingly relevant for higher 
latitudes) . 

We do not consider 7— ray and radio observations of the galactic center (GC) and galactic ridge since details of 
the DM profile are crucial for their interpretation. We also do not include 7— ray observations of dwarf spheroidal 
galaxies which are thought to have significant concentrations of dark matter. Purported evidence of this is that such 
satellite galaxies have survived several traverses of our galactic disk, stripping DM with each pass (which makes the 
use of conventional halo profiles difficult to justify). Interesting constraints have been obtained l39l. liol using HESS 
observations above 250 GeV of the Sagittarius Dwarf [4l[, and of the GC [42j and galactic ridge |43fTFor example, 
DM annihilation into dominantly t + t~ is disfavored as an explanation of the PAMELA and Fermi data for cusped 
halo profiles because of the accompanying large 7— ray flux from the decays of neutral pions [4fl |. However, DM decay 
into only t + t~ is still viable independently of the halo profile [40]. We mention this in anticipation of our results. 
We find that to be fully consistent with the Fermi and HESS e + + e~ flux while simultaneously explaining PAMELA 
data, a somewhat softer spectrum than obtained from DM annihilation or decay into e + e~ is required, as can be 
provided by the and t + t~ modes. 




FIG. 2: Best-fit x 2 values for DM annihilation into various charged lepton final states for several values of E B from a joint 
analysis of the PAMELA, Fermi 7— ray, Fermi e and HESS datasets which have 7, 18, 26 and 8 points, respectively. (e,fi,r) 
denotes a democratic final state with equal branching fractions into e + e~ , and t + t~ . The number of free parameters is 

8, and the number of degrees of freedom (dof) is 53, including two energy scale normalizations; see the appendix for details. The 
X 2 /dof is provided above each column. For E s = 150 GeV, all data other than the falling HESS spectrum are well reproduced 
for the fi + fi~ and t + t~ modes. An isothermal halo profile is assumed. 
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FIG. 3: Illustrative spectra for DM annihilation with E s = 150 GeV, 800 GeV, 1 TeV and 2 TeV into fx + fj,~ , e+e", (e, n, r), 
and t + t~, respectively. The DM-only contribution to 7— rays for the (J. + ft~ mode is tiny (because the boost factor is about 
25; see Eq. [7j) and lies outside the range of the figure. The corresponding \ 2 values can be found in Fig. [2] The background 
contribution for each case is displayed in the same line-type as for the signal. From top to bottom in the bottom-right panel, 
the electron background injection spectral indices are 70 = 2.38, 2.49, 2.57 and 2.63; see the appendix. The orange, light blue 
and dark green HESS data points (color coded according to the annihilation channel) are shifted within the energy calibration 
uncertainty with e = 1.3, 1.02 and 0.84, respectively. The HESS data corresponding to annihilation into (e, (X, r) have e = 0.98, 
and are not shown. For E s — 800 GeV, 1 TeV and 2 TeV, the Fermi e data shift downward (shown in black) with an energy 
calibration scale of 0.9, and for E s = 150 GeV (shown in gray), the corresponding value is 1.03. The shifts in the datasets 
are amplified by the E\ factor in the ordinate. The HESS and Fermi error bars have been expanded to approximately include 
systematic uncertainties (apart from the energy scale uncertainties). 



(i) Annihilating Dark Matter. The injection spectrum of high energy diffuse photons is 

^ = ^ BF ^ M ^; +C / -Ic(4^e;^)-e(^K SRF (4)^;^, (4) 

where the first term is the direct contribution from the annihilation process and the second term is the contribution 
from IC scattering of DM electrons (and positrons) with the ISRF at the same location. Here, <I> 7 is the photon flux 
per unit volume, BF is the boost factor assuming the canonical value (cry) = 3x 10~ 26 cm 3 s _1 , ^ 1 is the normalized 
photon spectrum from each annihilation event, 7i 7 / e is the interstellar photon and DM electron (and positron) densities, 
and aic(E' , E' e ;E^) is the differential IC scattering cross section for electrons of energy E' e scattering off photons in 
the ISRF with energy E^ and producing a photon with energy E y . A generic feature of annihilating DM is that the 
prompt photon flux depends sensitively on the DM halo distribution as p 2 . 
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(ii) Decaying Dark Matter. In this case, the flux of prompt photons, 



^ (5) 



T M DM dE 7 ' 

replaces the first term in Eq. ([4]) where T is the DM lifetime. Since the prompt photon flux depends linearly on the 
DM density, fewer photons originate from the GC than in the annihilation scenario. 

(in) Pulsars. Simulations indicate that the primary and 7— ray flux from pulsars have a luminosity ratio 
dN e ± / dE e ± : dN 7 /dE 7 between 0.2 and 0.5 [44[ for 7— rays detected by Fermi LAT. The primary 7— ray flux is a 
negligible component of the diffuse flux that is dominated by background contributions from IC scattering and 7r° 
decay. We include the contribution from IC scattering of the primary e flux emitted by pulsars. 

Results. We simulate the following energy spectra, and perform a likelihood analysis as described in the appendix: 

* _ for PAMELA 

£f ($ e + + $ e - ) for Fermi e ± and HESS (6) 
£^$ 7 for Fermi 7-ray 

where $ is the background plus signal differential flux of each species. 

In Fig. [2] we show the best-fit \ 2 values for DM annihilation into various charged lepton final states. The combined 
analysis includes the Fermi LAT, PAMELA, and HESS datasets. For PAMELA, we only analyze data above 10 GeV 
for which the solar modulation effects are negligible. We do not include WMAP haze data in our x 2 — analysis because 
it may be subject to unknown and potentially large systematic uncertainties. It is evident that the data prefer softer 
spectra. The Fermi e^ data and steepening of the HESS spectrum above 1 TeV are easily reproduced by heavy DM 
(Mdm — 2 TeV) annihilating to the t + t~ mode; see Fig. [3] The results for DM decay are only marginally worse, so 
we do not show them. 

Figure [4] shows a comparison of our three scenarios, all of which satisfactorily explain the experimental data. DM 
annihilation explains the data marginally better than DM decay. For DM annihilation, the p 2 enhancement near the 
GC yields a larger fraction of e ± that propagate a long distance and hence a slightly softer e ± spectrum that fits 
PAMELA and Fermi e ± data better. For DM annihilation and decay into we find the best-fit boost factor and 

lifetime to satisfy the approximate relations, 

BF = 43LE S - 38.9 , -J— = 2.29 + ^ , (7) 

10 Zb sec E s 

with E a in TeV. The exponential cut off in the pulsar electron spectrum can reproduce the HESS data very well. 
Notice that the the quality of the fit is only weakly dependent on E s for pulsars. 
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FIG. 4: A comparison of the best-fits for DM annihilation, DM decay and pulsars. In the DM scenarios, the fj> + mode is 
marginally preferred over the r + r~ mode, except for E s = 2 TeV. 
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FIG. 5: Illustrative spectra for DM annihilation into fj, + fj,~ (E s = 1 TeV), DM decay into t + t~ (E s = 2 TeV), and pulsars 
with E s = 1 TeV. Note that the e spectra for DM annihilation into r + r~ with E s = 2 TeV are identical to the DM decay 
case. The pulsar-only contribution to 7— rays lies slightly outside the range of the figure. The corresponding \ 2 values can be 
found in Fig. [4] For the DM decay case, the green HESS data points are systematically shifted to lower energies with e = 0.84. 
For pulsars, the light blue HESS data points have e = 1.0. The HESS data corresponding to DM annihilation have e = 0.8, 
and are not shown. For the DM cases, the Fermi e ± data shift (shown in black) with an energy calibration scale of 0.9, and for 
pulsars (shown in gray), the corresponding value is 1.0. The electron backgrounds have 70 — 2.61. 



In Fig. O we show representative best-fit spectra from the joint analyses of the Fermi, PAMELA, and HESS data. 
The mid-latitude preliminary Fermi 7— ray data reside in the energy range where the galactic background significantly 
dominates the DM contribution. Above about 50 GeV, the DM signals dominate, and may be resolvable with future 
Fermi LAT data up to 300 GeV. Even with future data, the isotropic EG component is significant at mid-latitudes 
and needs to be carefully investigated before robust conclusions can be drawn. From Fig. [5] we see that the 7— ray 
contribution from pulsars is negligible. 

In Fig. [5J we show how DM annihilation signals depend on the DM halo profile. We single out the annihilation 
case because of its stronger dependence on p. For a cuspier profile, a larger fraction of e are injected near the GC so 
that more synchrotron emission occurs as they propagate through the galaxy. This is evident in the upper right-hand 
panel of Fig. [5J We see that there is a marginal preference for cusped profiles essentially because the combination of 
PAMELA and Fermi e ± data require a not-too-hard spectrum. The different profiles produce similar spectra, 
and require boost factors of about 400 with a difference of less than 3% among them. 

Although Fermi LAT does not corroborate the 7— ray excess in EGRET data at mid-latitudes, it is still quite 
possible that it may observe an excess at the GC; see Fig. [7] In this region of the sky the DM signal dominates the 
EG background. In contrast with the isothermal profile, the cusped profiles enhance the flux in this energy range, 
signaling DM annihilation and a cusped halo profile. 
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FIG. 6: Impact of different DM halo distributions for DM annihilation into with E s = I TeV. The cuspy Moore, Einasto 

and NFW profiles yield higher contributions to the diffuse 7 and synchrotron spectra than the isothermal profile. The e ± 
spectra for the different profiles overlap and have similar boost factors. The HESS and Fermi energy calibration scales are 
e = 0.8 and 0.9, respectively, and are the same for all profiles. The WMAP haze data from the outer [22J and inner [23] galactic 
regions are shown at a representative frequency and only with statistical uncertainties, are not included in the \ 2 ~ analyses, 
and are unexplained by our scenarios. The ATIC data are not included in the fit, but are displayed for comparison. The 
corresponding x 2 values show that the data do not favor a particular profile. 
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FIG. 7: Diffuse gamma-ray spectra in the GC region \l\, \b\ < 5°. The left-hand panel shows spectra corresponding to the DM 
annihilation, decay and pulsar cases in Fig. [5] where DM assumes the isothermal halo profile. The right-hand panel shows the 
spectra corresponding to the DM annihilation cases with different profiles in Fig. [6] EGRET data [191 ] are plotted to illustrate 
that a significant diffuse gamma-ray excess can be expected from cuspy profiles in the GC region although preliminary Fermi 
LAT results see no excess in the mid-latitude range. 



Summary. Three popular explanations have been proposed for the sharp rise in the positron fraction in PAMELA 
data: DM annihilation, DM decay and pulsar sources. Associated with these possibilities are 7— ray signals. We 
deliberately avoided existing observations of the galactic center since they are sensitive to the details of the halo 
profile, and focused on the diffuse 7— ray background from latitudes above 10°. Preliminary results from the Fermi 
LAT do not confirm the diffuse 7— ray excess at mid-latitudes reported by the EGRET observatory. With these 
data that are consistent with the expected background, we entertained the possibility that final state radiation and 
inverse Compton scattering (which are the primary sources of 7— rays) may help constrain DM as an explanation 
of the PAMELA, Fermi e ± and HESS data since these data require copious production of positrons and electrons. 
However, we found that this is not the case for DM annihilation/decay into charged leptons; our choice of leptophilic 
annihilation/decay was dictated by the nonobservance of an excess in the p/p flux ratio measured by PAMELA. We 
emphasize that the preliminary Fermi results are in a 7 energy range (< 10 GeV) where the galactic and extragalactic 
backgrounds overwhelm the DM contribution. With upcoming observations up to 300 GeV, Fermi data promises to 
be far more constraining. Dark matter with E s between 800 GeV and 2 TeV is easily consistent with all 4 datasets 
for both annihilation and decay. The data prefer the softer spectra that arise from the /i + /i~ and t + t~ modes; see 
Fig. 0J The boost factors required for DM annihilation and lifetimes for DM decay follow Eq. ([7]) for the mode. 
Overall, pulsars have the least tension with the data, primarily because the HESS data show a gradual fall-off in 
the energy spectrum, which is characteristic of pulsar sources. To draw these conclusions we varied the backgrounds 
within reason, accounted for available systematic uncertainties, and assumed an isothermal profile, ft is no surprise 
that with an isothermal profile, the WMAP haze could not be accommodated. 

For DM annihilation, the source term depends on the DM density as p 2 , so effects of the halo distribution are 
especially relevant. We considered 4 different profiles and found that the 4 datasets when analyzed jointly show no 
appreciable preference for the degree of cuspiness. The observance of a significant gamma-ray excess from the galactic 
center but no excess at mid-latitudes is a striking signal of DM annihilation and a cusped DM halo. 

Acknowledgments. We thank M. Medvedev for a discussion. This work was supported by DOE Grant Nos. 
DE-FG02-04ER41308, DEFG02-95ER40896, DE-FG02-84ER40173 and DE-AC02-06CH11357, by NSF Grant No. 
PHY-0544278, and by the Wisconsin Alumni Research Foundation. 

Appendix. Numerical procedure. To fully simulate the e ± and photon spectra, we use the GALPROP package [28| 
with the ISRF profile [45]. The DM 7— ray and e ± injection spectra are calculated using DMF1T and Mi- 
crOMEGAs [13], respectively. We run GALPROP in the 2-dimensional mode and use the "conventional" (ID:500180) 
model [48| for the galactic set-up and cosmic ray source/injection parameters; however, we vary the injection spectral 
index 70 for the electron background (parameterized as a power-law E~ 10 ) in the range 2.2 — 2.9. This model produces 
a p/p spectrum that agrees with PAMELA data above 10 GeV. The cylindrical diffusion zone has a radius of 20 kpc 
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and a height of 4 kpc; the magnetic field that determines the synchrotron radiation intensity is given by 



B = B Q e ~ e" , (8) 

with a default local strength Bq = 5 pG and scales r s = 10 kpc, z s — 2 kpc. The final 7 spectra are calculated from 
the GALPROP created sky-maps by averaging the IC, FSR, bremsstrahlung, and ir° decay contributions over specific 
ranges of galactic coordinates. The galactic and source spectra are generated separately and then combined with a 
varying normalization factor for the latter. 

Cosmic ray propagation is governed by the energy loss - diffusion equation, 

^-D(E)- V H-d E (D p (E)-$) = Q, (9) 

where Q is the injection source term for each species. The energy loss coefficient D p (E) is dynamically evaluated 
and for e ± it is mainly determined by bremsstrahlung and IC scattering. The spatial diffusion parameter D(E) is 
parametrized {49i| as 



s 

2 r 



where j3 = v/c = 1 for cosmic rays of interest. R(E) and Rq are the particle rigidity at energy E and the reference 
rigidity, respectively. For e^, R(E) = E. To account for the uncertainty in cosmic ray diffusion, we vary the 
diffusion parameters (Do, Rq, 5). Since the various D(l GeV) fall in the range (3 — 5) x 10 28 cm 2 /s to fit cosmic ray 
nuclei data [4!|, we impose this as a constraint in our analysis. We restrict S to be above 1/3, the value expected 
for a turbulent magnetic field with a Kolmogorov spectrum. For each E Sl we also vary the total background e~ 
normalization, and the BF (for DM annihilation), or T (for DM decay), or N (in Eq. [2] for pulsars). 
For the DM halo, we considered the following halo distributions: 



p{r) 



P& 1 t+(flr s yi Isothermal, r s = 5 kpc [29] 

Po " 



£ (l + £) " NFW, po = 0.26 GeV/cm 3 , r s = 20 kpc [30] 

p©exp{-f [(r a - rg)/r«]} Einasto, a = 1.7, r s = 25 kpc [31] 



k Pe (^) 1 - 16 (t^t) 1 ' 84 Moore, r s = 30 kpc [32] 



where the local DM density is pq — 0.3 GeV/cm 3 for all profiles, and r s is the size of the DM halo. The spatial source 
of e , 7 is obtained by convolving the halo distribution with the injection spectra. 

Likelihood analysis. We consider Fermi, PAMELA and HESS as independent sets and sum the x 2 contribution from 
each set. In obvious notation. 



/ fth fex\2 

■2 X < X ' Ui h_J_ , 12 x 

A (Af ex ) 2 

experiments i 



where / are the spectra in Eq. ([5]). Since bin sizes are available for PAMELA and Fermi e ± data, we average over 
each bin to obtain the theoretical expectation f th . 

The HESS data have a 15% energy calibration uncertainty. To account for it we fit an additional parameter 
e = 1 + that scales a data point from (E, E 3 ^) to (eE, e 2 E 3 ^), while the number count E^ remains fixed. 
With f(E)=E 3 ^, 

2 M _ V (/ tft (^)-e 2 /f^ S ) 2 . (1-e) 2 m] 
XhessI 6 ) ~ 2^ (e 2 Af l HESS ) 2 + (A e )2 ' [i6 > 

where Ae = 0.15 and Aff 1ESS includes the additional systematic uncertainty (from modeling hadronic interactions 
and the atmosphere) in quadrature with the statistical uncertainty [T3|. The Fermi e ± data have a ^^o^ energy 
calibration uncertainty which we treat in analogy with the HESS calibration uncertainty, by including an additional 
energy scale parameter, except for one important difference: We allow this parameter to vary only within the la range 
of the calibration uncertainty so that the Fermi e data do not rigidly shift too much and become discordant with 
data from AMS, HEAT, ATIC, BETS and PPB-BETS which are consistent with each other between 10 GeV and 
100 GeV. We have checked that the Fermi dataset shifts significantly unless low energy data from other experiments 
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are included in the analysis; the low energy data from ATIC, for example, anchor the Fermi energy calibration so that 
the scale parameter does not deviate too far from unity. 
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